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We provide a numerical study showing that a bottom reflector is indispensable to 
achieve unidirectional emission from a photonic-crystal (PhC) nanolaser. First, we 
study a PhC slab nanocavity suspended over a flat mirror formed by a dielectric 
or metal substrate. We flnd that the laser's vertical emission can be enhanced by 
more than a factor of six compared with the device in the absence of the mirror. 
Then, we study the situation where the PhC nanocavity is in contact with a flat 
metal surface. The underlying metal substrate may serve as both an electrical 
current pathway and a heat sink, which would help achieve continuous-wave lasing 
operation at room-temperature. The design of the laser emitting at 1.3 /xm reveals 
that relatively high cavity Q of over 1,000 is achievable assuming room-temperature 
gold as a substrate. Furthermore, linearly-polarized unidirectional vertical emission 
with the radiation efficiency over 50% can be achieved. Finally, we discuss how this 
hybrid design relates to various plasmonic cavities and propose a useful quantitative 
measure of the degree of the 'plasmonic' character in a general metallic nanocavity. 
@ 2011 Optical Society of America 

OCIS codes: 140.5960, 140.3945, 230.5298, 240.6680. 
1. Introduction 

Spontaneous emission of a dipole emitter can be altered by the presence of a metallic or dielectric re- 
flector. For instance, the spontaneous emission may be completely inhibited inside an appropriately 
designed optical cavity. Purcell provided the flrst quantitative analysis on the emission dynamics 
of an atomic dipole placed inside a cavity characterized by its quality factor {Q) and mode volume 
(V^). [1] It has been now well established that the effect of the spontaneous emission modification 
will be more pronounced with the higher Q/V ratio. [2] 

On the other hand, researchers in the field of semiconductor lasers have searched various ways 
to achieve the so-called 'thresholdless' laser in the context of spontaneous emission control by some 
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form of high QjV cavities. [3] The concept of photonic-crystal (PhC) [4,5] has revolutionized the 
development of such high QjV lasers, enabling further miniaturization in device size and reduction 
in threshold power. The periodic arrangement of dielectrics can result in a forbidden frequency 
region within which any electromagnetic mode cannot propagate, a property now called the 'pho- 
tonic band gap (PBG).' [6] In fact, the presence of such an energy band gap in a one-dimensional 
(1-D) periodic structure (Distributed Bragg reflector, DBR) [7] has been known for some time and 
already applied to the design of the vertical-cavity surface-emitting laser (VCSEL). [8] The artificial 
material possessing PBG enables us to confine electromagnetic energy in a volume smaller than the 
associated wavelength of light. The first nanocavity laser was achieved based on a thin dielectric 
membrane with periodically arranged air-holes in two-dimensions. [9] Laser gain was provided by 
multiple layers of quantum wells embedded in the middle of the slab. Though Q of the initial laser 
cavity was below 500, much progress has been made toward higher Q and smaller V and this has 
been a major research topic until very recently. [10-12] However, the importance of out-coupling 
efficiency in these wavelength-scale emitters has for some time been neglected. 

It is interesting to notice that the basic architecture of PhC nanocavity design has not changed 
much since the original air-suspended thin-slab geometry. [13] The use of such a thin slab appears 
to be indispensable to maximize the size of the in-plane PBG. [14] For a high refractive index 
semiconductor slab, a typical thickness of the PhC slab is 0.5 a, where a is the lattice constant. 
Yet in reality, most PhC slab nanocavities are suspended over a dielectric (or metallic) substrate, 
as depicted in Fig. 1(c). Thus, they are no longer isolated but experience feedback from their 
environment. It has been pointed out that even cavities with Q in the range of 50,000 can 'see' the 
underlying mirror surface, leading to severe modification in the far-field radiation profile. [15, 16] 
This phenomenon, as is to be discussed in the following section, has a strong analogy to the 
aforementioned cavity quantum electrodynamics (QED) example of the point dipole source near a 
plane mirror. [17] We will show that the combination of a highly reflective bottom mirror and a 
proper gap size between the PhC slab and the mirror can lead to enhanced far-field directionality 
from the PhC cavity by more than a factor of 6 in comparison with the PhC nanocavity in the 
absence of the bottom reflector. This result should be of practical importance for the various PhC 
based light emitters including nanolasers and single-photon sources. [16, 18, 19] 

Recently, the realization of an electrically-pumped nanolaser has drawn renewed attention of 
many researchers. The first current injection PhC nanolaser was demonstrated in 2004, in which a 
submicron-sized dielectric post was introduced right below the laser cavity as a means to deliver 
and confine electrical current. [20,21] However, in order to obtain reasonably high Q for lasing, the 
post structure needs to be made very thin (diameter < 500 nm) and long (typically l^m). This 
could lead to unusually high electrical resistance of over IkQ, which has been a major bottleneck to 
achieve continuous- wave (CW) operation of a laser at room-temperature (RT). To circumvent this 
issues, laterally-doped p-i-n structures were proposed recently. [22-24] CW lasing has been achieved 
by external cooling at an ambient temperature of ~ 150 K. [23] Note that optically-pumped RT- 
CW lasing operation of a PhC nanolaser was already demonstrated by Nozaki, et al, [25] however, 
realization of its electrical counterpart has been a severe challenge. 
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Fig. 1. Evolution of photonic crystal nanolaser: From vertical-cavity surface-emitting 
laser to hybrid metal-photonic crystal laser. 

Therefore, we believe now is the time to reconsider the PhC slab design itself. The concept of 
using a bottom mirror can naturally lead to the following question - What would happen if we make 
the air-gap size zero by placing the whole PhC slab directly on the mirror? We show that a flat 
metal substrate (See Fig. 1 (d)) can be used effectively to achieve good vertical confinement of the 
cavity mode. [26] Furthermore, metals are very good conductors for both electrical current and heat 
dissipation. This may imply that the aforementioned difficulties in building current-injection PhC 
nanolasers could be mitigated using this hybrid metal-PhC slab design. We shall provide numerical 
simulation results on the metal-bonded PhC nanocavities in Sec. 3. Since the PhC structure is now 
in contact with the metal substrate, we can expect certain 'plasmonic' effects in this hybrid design. 
In fact, it may not be easy to distinguish between 'photonic' and 'plasmonic' when both characters 
coexist in the same structure. Therefore, we will discuss the relationship of this newly proposed 
design to various 'plasmonic' cavities. We will propose a quantitative measure of the degree of 
the 'plasmonic' character, ^ plasmonicity\ which would provide a useful guideline in the design of 
'plasmonic' cavities. 

2. A PhC nanocavity nearby a bottom reflector 

2. A. Cavity QED analogy 

Before going into the details of the hybrid metal-PhC nanocavity, let us first consider a PhC 
nanocavity suspended over a flat mirror as depicted in Fig. 2(a). The situation is directly analogous 
to a point dipole emitter near a plane mirror, which is a well-known cavity QED example. [17] 
Certainly, there are differences between the PhC nanocavity and the dipole emitter. The PhC 
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Fig. 2. (a) A photonic crystal nanocavity is suspended in air above a flat mirror (a 
bottom reflector) . The radiative decay rate (7) of the nanocavity mode can be tuned 
as a function of the air-gap size, (b) The design of the photonic crystal nanocavity. 
Here, two air-holes facing each other are enlarged by Rp = 0.05a. Other parameters 
are as follows: the slab thickness (T) = 0.9a, the modified hole radius {Rm) = 
0.25a, and the background hole radius {R) = 0.25a. The lattice constant of the 
photonic crystal is denoted as 'a' throughout this paper, (c) Electric-field intensity 
distribution (|Ep) of the deformed hexapole mode detected in the middle of the 
slab (z = 0). (d) Normalized decay rates (7/70) of the deformed hexapole mode 
as a function of the air-gap size. Perfect electric conductor (PEC), gold, and a 
dielectric of the same refractive index of 3.4 as the slab material are considered as 
a bottom reflector. In the case of a gold mirror, we assume emission wavelength to 
be ~ 1.3;um with a = 450 nm. Drude model parameters are as follows; eoo = 10.48, 
ujp = 1.38 X lO^^rad/s, and 7™ = 1.18 x lO^^rad/s. (e) The resonance frequency also 
changes as we vary the air-gap size. 
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nanocavity is much bigger in size than the atomic dipole emitter. The radiation pattern from this 
hypothetical emitter does not resemble to that of the simple dipole emitter. In general, the complete 
description requires higher-order multipoles including both electric- and magnetic-multipoles. [27] 
Thus, we can view this generalized light emitter as a sum of point-like multipole emitters. Similar 
modifications in the decay rate and the radiation pattern from the PhC nanocavity are expected. 
In the absence of a bottom mirror, we find that Qq of the deformed hexapole mode (Fig. 2(b)) 
is ~ 15,000 (subscript for the cavity with no bottom reflector). [15] To make our analysis more 
analogous to [17], we translate Qq into the radiative decay rate using the relation 70 = uj/Qq. 

We have performed 3-D finite-difference time-domain (FDTD) simulations to study how 7 changes 
in the presence of a bottom reflector. In Fig. 2(d), we plot normalized decay rate (7/70) of the 
deformed hexapole mode as a function of the air-gap size (d). Three different types of the bottom 
reflector are considered in this study - perfect electric conductor (PEC), gold, and dielectrics. 

First, let us focus on the ideal mirror case (PEC). As expected, we can observe modulations 
of 7 as a function of d. Both enhanced decay (7/70 > 1) and suppressed decay (7/70 < 1) can 
occur depending on a specific d value. Interestingly, we have found that the modulating features 
appear to have a certain periodicity. For example, the two consecutive peak positions (one at 1.75 
a and the other at 3.5 a) are separated by ~ 0.5A. Here, it should be noted that, la is equal to 
~ 0.29A for the cavity displaced sufficiently far from the underlying mirror, where A is the emission 
wavelength measured in vacuum. The observed ~ 0.5A periodicity reminds us of the two-beam 
interference condition, where the bottom reflector can reverse the downward propagating beam to 
the upward direction to produce the interference. We will develop a more rigorous model in the 
following subsection 2.B. 

Second, the case of a gold reflector is considered. We use the dielectric function of gold at RT for 
emission wavelength of ~ 1.3 ^m.. [28] Slightly increased decay rates in comparison with the previous 
ideal mirror case are partly due to the additional absorption and the slightly lower reflectivity of 
gold. The absorption effect becomes more severe for smaller d. Indeed, the decay rate becomes 
noticeably different from the PEC case when d < 0.5a. 

Third, when we replace the gold mirror with a simple dielectric substrate with refractive index 
of 3.4, the periodic modulation of 7 becomes much weaker than the previous two cases. This is not 
surprising because the mirror reflectivity is now only about 30%. For d < 1.2a, we observe much 
more enhanced decay rate — at d = 0.5a, 7/70 is about 38. Such a dramatically enhanced decay 
rate is mostly due to the enhanced tunneling loss through the bottom substrate and the TE-TM 
coupling loss in the horizontal direction. [15] In fact, light confinement mechanism in the in-plane 
directions is no longer perfect as we break the vertical symmetry of the PhC slab. [14] 

In all three cases, we observe the break-down of the PEG as d — t- 0. However, it is much more 
severe in the case of a dielectric mirror. We would like to note that 7/70 converges to 1 in the opposite 
limit of d — 7- 00, though it has not been shown explicitly in those plots. Finally, in Fig. 2(e), it 
can be seen that the cavity resonance blue-shifts (Acjbiue) as d decreases. Such energy-level shifts 
are also observed from the aforementioned cavity QED example. [17] One may develop a similar 
perturbative approach based on electric- and magnetic-multipoles to explain the observed Awbiue- 
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However, arguments based on the electromagnetic variational theorem [6] would suffice to explain 
differences in At^biuc in the three different cases — the dielectric mirror case shows the smallest 
Acjbiuc due to the more efficient overlap of the electric-field energy with the dielectric mirror region. 

2.B. Enhancing energy directionality: Planewave interference model 

As mentioned previously, the interference of electromagnetic waves is mediated by the bottom 
mirror to produce the observed decay rate modulation. Remember that the optical loss of the 
nanocavity is closely related to its far-field radiation pattern, therefore, we expect that the far- 
field radiation pattern should undergo similar modifications. By using 3-D FDTD, we can directly 
simulate a far-field radiation pattern, dP{6,(j))/dO,, of the PhC nanocavity, which represents emit- 
ted power (dP) within a unit solid angle (dO). There exists corresponding wavevector component 
k for each angular direction {6,(j)). Thus, any radiation pattern can be decomposed in terms of 
planewaves ^i^.((iP/dr2)(k). This planewave decomposition provides an alternative way to the mul- 
tipole expansion method. From now on, we will focus on planewaves with kx = and ky = 0, since 
we are interested in the vertical directionality of laser emission. 

Figure 3(b) describes how the complex 2-D slab nanocavity is simplified in the spirit of the 
planewave decomposition. The perforated PhC membrane is approximated as a uniform dielectric 
slab with an effective refractive index. [29] We will deal with the PEC mirror case only, where any 
wave incident upon it will undergo a vr-phase shift. Detailed calculations have been described in our 
earlier publication, [15] here we would like to summarize the essential ideas. First, the final result 
will depend only on the effective index of the slab, nefj, the phase thickness of the slab, (p, and the 
phase thickness of the air-gap, tp. Once is determined, we can derive coefficients of amplitude 
reflection and transmission for a single dielectric interface such that rg = (ngfj — l)/(nefj -|- 1) and 
to = (1 — ^o)^^'^- Second, all of the scalar waves will be treated as complex numbers in the form 
exp{ikzZ — iojt). Any wave traveling across a phase thickness of (f)Q will gain the same amount of 
phase. Third, we assume that there are two 'seed' waves, one on the top of the slab and one at the 
bottom of the slab. Both waves are assumed to propagate in the opposite directions with the equal 
amplitude and phase. This assumption can be justifled by the fact that the original PhC nanocavity 
mode is symmetric with respect to the z = plane. Fourth, the wave initially propagating in the 
downward direction will be redirected upward by the bottom mirror. In Fig. 3(b) S denotes the sum 
of all such waves finally detected in the far-zone at = 0. During this complex process, a fraction 
of the energy can couple to the PhC nanocavity mode. If this coupling occurs, then the resonant 
mode can produce wavevector components with kx 7^ or ky 7^ 0. In this model, we assume that 
this coupling process is negligible and the in-plane wavevector components are conserved. This 
assumption is analogous to the weak-coupling regime in the cavity QED [17] — once a photon 
leaves the emitter (the PhC resonant mode), it does not strongly interact with the original emitter. 
The photon will 'see' only mirror boundaries of an external cavity (the bottom mirror and the 
uniform dielectric slab) whose Q is much smaller than that of the emitter. It should be noted that 
similar approaches have been applied to the calculation of spontaneous emission inhibition and 
enhancement of a dipole emitter in a cavity. [30] 
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Air-gap size (9) 



Fig. 3. (a) Vertical emission enhancement factor (W) obtained by the planewave 
interference model. Both the air-gap size and the slab thickness are varied and the 
results of varying slab thicknesses were shown as multiple curves as a function of 
the air-gap. We have assumed the effective refractive index of the slab (^^e//) to 
be 2.6, which will result in tq = {n(,fj — l)/(ne// -|- 1) ~ 0.44. (b) A schematic of 
the model for the photonic crystal nanocavity suspending over a bottom reflector. 
The perforated slab is replaced with a uniform dielectric slab with n^fj and the 
underlying mirror is assumed to be PEC. 
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With all the assumptions above, the relative vertical enhancement can be written as |1 + Sp, 
which is normalized to the result in the absence of the bottom reflector. Here is the final result of 
the summation. 

We plot as a function of the air-gap size {(p), where we assume n^g to be 2.6 in consideration of 
the effective surface coverage ratio of ~ 1 - [2tt R"^ / (VSa^)] for the PhC slab with R = 0.35a. Fig. 3 
shows the results for various (p values (the slab thickness). When (p satisfies the 'slab resonance' 
condition (0 = m x vr, where m is an integer), the slab will be 'transparent' for any wave incident 
upon it. Under this condition, the redirected wave from the bottom mirror does not 'see' the slab, 
which result in the simple two-beam interference. W modulates between and 4 by varying the 
air-gap size. The effective thickness of a PhC slab is typically chosen to be around T ~ A/(2ncff)' 
in order to maximize the PBG, [14] which happens to be near the slab resonance condition. This 
implies that near 100% transmission is possible by slightly tuning the emission wavelength toward 
the exact slab resonance condition. [29] Then, the subsequent optimization of the air-gap size 
can bring up the vertical emission intensity by up to a factor of 4. Here, it is worth mentioning 
that the enhancement factor is a quantity proportional to the decay rate. Thus, the fact that 
Wmax = 4 > 1 -|- 1 does not violate the conservation of energy — most of the electromagnetic 
energy, U{t), is stored in the cavity and it will experience exponential decay in time such that 
U{t) = U{to) exp[7(t — to)], where the decay rate, 7, is proportional to W. Conversely, if the air-gap 
size is chosen to be I.Stt, then, for an ideal 1-D system, we can completely quench the radiation 
to the vertical direction, resulting in infinitely large Q. Similar unbounded Q behavior has been 
found from the second F-point band-edge state in a 2-D honeycomb-lattice PhC slab [31] and the 
F-point band-edge state in a 2-D triangular-lattice PhC slab, [32] with and without using a bottom 
reflector, respectively. 

Now let us turn to our original problem and find out the condition that maximizes W. Surpris- 
ingly, W can increase up to about 9 by tuning the slab thickness near the 'slab antiresonance' 
{(j) = (m + 0.5) X vr) as shown in Fig. 3(a). [15] This phenomenon can be understood by the fact 
that the slab becomes highly reflective (reflectance can be as high as 70% due to the absence of 
the slab resonance). [29] Thus, the slab and the bottom mirror system constitutes a certain vertical 
resonator, which enhances the density of photon modes near = ky = point. In the follow- 
ing subsection, through rigorous 3-D FDTD simulations, we will show that over six-fold vertical 
emission enhancement can be achieved in the deformed hexapole mode example. 

2.C. Enhancing energy directionality: FDTD 

The planewave interference model predicts that vertical emission can be enhanced by more than a 
factor of 4. Here, we would like to test the validity of the model by using rigorous FDTD simulation. 
In our previous paper, we have adopted the near- to far-field transformation algorithm for more 
efficient and fast simulation. [15] However, for more accurate results, we will perform direct FDTD 
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Fig. 4. FDTD simulated far-field emission profiles from the deformed hexapole mode 
shown in Fig. 2. Far-field patterns detected over the hemispherical surface are trans- 
formed into the 2-D plane by using a simple mapping defined hy x = 9 cos (p and 
y = 9 sin (p. Numbers represent the air-gap size normalized to the emission wave- 
length of the reference cavity (cxo) in the absence of the mirror. 
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simulations in a very large computational domain with {L^ x Ly x L^) > (6A x 6A x 6A). We 
take the same deformed hexapole mode shown in Fig 2 as an example and assume a PEC bottom 
mirror. Note that we will vary only the air-gap size while we keep the slab thickness constant at 
0.9a. First, far-field emission profiles, (^), are obtained by detecting the radial component of the 
Poynting vectors over a hemispherical surface whose radius is larger than 3A (See Fig. 4). These 
far- field patterns are then normalized by the original far- field pattern, fQ(6,(p), obtained in the 
absence of the bottom reflector (See Fig. 4). By varying the gap size, it is found that the internal 
electromagnetic energy, U{t), does not show noticeable change within the FDTD time needed for 
the far-field simulation. Thus, errors associated with small variations in U (t) will be less than 
1%. We present normalized far-field patterns, f{0,4>) = f / fo, in Fig. 5(a), where we use a simple 
mapping defined hy x = 6cos(j) and y = 6sin(j). Thus, the center corresponds to 6 = and the 
angle 6 will be proportional to the radial distance from the origin. Clearly, we can have more than 
six- fold enhancement by choosing the air-gap to be near lA. 

In Fig. 5(b), we plot f{0 = 0) as a function of air-gap size. The three solid lines are from 
the previous planewave interference model, for (p = l.lvr, 1.27r, and l.Svr, respectively. Generally, 
the curve for (j) = l.Svr shows good agreement with the FDTD results. Interestingly, the model 
also predicts the asymmetric behavior around l.OA as is confirmed by the FDTD result — the 
enhancement value decreases more slowly on the left side than on the right side. Thus, both the 
FDTD result and the planewave interference model confirm that the vertical emission can be 
enhanced by more than a factor of 6 by employing a highly reflective mirror at the bottom. 

2. D. Applications 

These findings will be of significant importance to the design of various PhC based light emitters, 
such as nanolasers and single-photon sources. [16, 18, 19] The emission wavelength is a rather fixed 
property by the active material embedded in the slab. In most wafer designs, the final selective 
wet-etching process will leave a flat surface of a dielectric substrate, which will serve as a good 
bottom mirror as we have shown in Fig. 2. Furthermore, in GaAs material systems, we may include 
a GaAs/AlAs DBR whose resultant reflectivity will be over 98%. To optimize the vertical direc- 
tionality, the following two criteria should be considered for the wafer design. 1) The thickness of 
the sacrificial layer (which is to become the air-gap) should be designed to be equal to the emission 
wavelength. 2) The PhC slab thickness should be chosen to be near the slab antiresonance. How- 
ever, in consideration of the PEG size, the slab thickness may be limited by 1.0a ((/)max ^ l.Svr). It 
should be noted that, the enhanced directional emission has been verified from InP/InGaAsP PhC 
nanolasers by direct measurement of far- field emission profiles. [SS] 

3. A PhC nanocavity on a metal substrate 

As mentioned previously, the planewave interference model may not be applicable to the case where 
the air-gap size is below A/2. However, we have seen already what would happen in the limit of 
zero air gap from Fig. 2(d,e). In the case of the dielectric mirror, we cannot find any mode whose 
Q is higher than SO. However, a relatively high-Q (> 1,000) mode is found by assuming the PEC 
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mirror. The case of the gold mirror is not as good as the case of the PEC mirror — Q ^ 150 has 
been obtained. Though this Q value is already comparable with typical Q values from previously 
reported metallic nanocavities, [34-38] we will propose a method to bring it up to several thousands. 
In this section, we will explore another opportunity with this new metal-PhC nanocavity design 
(Fig. 6(a)) for a practical nanolaser. The fact that a metal substrate is in direct contact with 
the PhC slab structure may mitigate aforementioned difficulties in building electrically-pumped 
nanolasers, namely, excess electrical and thermal resistance. Moreover, the metal may work as a 
good bottom reflector so that the vertical directional emission could be enhanced. This hybrid 
metal-PhC design will also be an important building block in the field of plasmonics in the context 
of building efficient plasmonic lasers. [39] 

3. A. Q, V, and Purcell factor Fp 

Before beginning with our discussions on the result in Fig. 6, let us clarify the definitions of Q, V, 
and Fp and several energy related quantities, especially in the context where dispersive metals are 
involved. 

Qtot is defined through the decay rate of the total electromagnetic energy contained in the cavity, 
UEiviit), such that 



UEM{t) = Uem{0) exp 



UJ 

-t 



(2) 



Qtot 

Here, UEM{t) is the sum of the electric-field energy, UE{t), and the magnetic-field energy, UM{t) 
stored in the PhC nanocavity mode. UE{t) and UE{t) can be defined by the energy density functions 
UE{T^,t) and itA,/(r,t), respectively. 



UE{t) = / d-'ruE{r,t) (3) 
Jv 

Unit) = I d^ruM{r,t) (4) 



V 

In the case where dispersive media are involved in the energy calculations, special care must be 
taken. In this paper, we will use the following definitions for ue and um^ as has been noted by 
Chang and Chuang. [40, 41] 

' d{ijje) 



UE{r,t) = -yRe 



du 



(E(r,t)-E(r,t))^ (5) 



UM{r,t) = |Re[e(L^)](E(r,t)-E(r,t))y (6) 

Here, the bracket (• • ■)rp denotes time average over one optical period. If r lies in a normal dielectric 
medium of e^, then the above expressions will take the familiar forms of n£;(r,t) = UM{T^,t) = 
(eoed/2) (E(r, t) • E(r,t))^, where the equality between the electric-field energy and the magnetic- 
field energy comes from the fact that the resonant mode is harmonically oscillating in time. [42] 
If r lies in a Drude medium where em{^) = ^oo — ^p/{^'^ + i'Jm^), we obtain 




UE{r,t) = 6oc+ ^2 , 2 2 (E(r,t)-E(r,t))r (7) 
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(a) Normalized far-field pattern 
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Fig. 5. (a) The far- field emission profiles shown in Fig. 4 are normalized by the 
reference far-field pattern (oo in Fig. 4), where white regions denote values > 6.33. 
(b) We extract ^ = components from the normalized far-field patterns and plot 
them together with the theoretical curves obtained by the planewave interference 
model. 
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(E(r,t) •E(r,t))^ 
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Throughout this paper, we adopt the cavity QED definition for the mode volume V, which, using 
the above definition of Uem, can be written as 



where a single emitter is assumed to be located in a medium of a refractive index n. 

Now consider the dipole mode in the PhC cavity in contact with the gold substrate as shown in 
Fig. 6(a,b). As mentioned previously, Qtot cannot be larger than 200 in the case of the deformed 
hexapole mode in a thin slab (slab thickness T = 0.9a). Traditionally, the PhC slab thickness is not 
chosen to be much larger than 1.0a, because the in-plane PBG begins to close at this value. However, 
Tandaechanurat, et al. have shown recently that relatively high Q can be obtained from the PhC 
dipole mode suspended in air after the PBG closure. [43] They have found that Qtot can increase 
up to ~ 10, 000 at T ~ 1.4a. In Fig. 6(c), we plot Qtot and y as a function of the slab thickness. We 
have increased the slab thickness up to 6a ~ 1, 800 nm. Surprisingly, Qtot shows a rather monotonic 
increase and a value of 3,000 can be obtained when T > 1, 500 nm. Qtot niay be further improved 
by inserting a thin low refractive index layer between the PhC slab and the bottom mirror, [44] 
which will reduce the optical overlap with the underlying gold. What is interesting here is that 
Qtot seems to increase indefinitely with increasing T, although the in-plane PBG is completely 
closed. This unusual behavior may be understood based on the waveguide dispersion along the z 
direction {oo-kz diagram) — the structure can be viewed as a PhC fiber with a finite length. [6,45] 
In the case of a thin slab, is not a well-defined quantum number. However, as T increases, k^ 
of the fundamental mode can be more and more precisely defined in accordance to the uncertainty 
relation between Az and /\kz- In fact, kz of the fundamental slab mode will converge to zero in 
the limit of infinite T. This fcz = point corresponds to the ideal 2-D limit in the uo-kz diagram. 
Thus, the mode can be more and more confined as T increases. It should be noted that, with a 
certain large T, radiation loss will occur mostly in the horizontal directions {x-y plane) rather than 
in the vertical direction {z). [46] This is due to the presence of the zero group velocity dispersion 
at fc^ = point. The slow group velocity mode will effectively reduce the scattering loss at the top 
surface of the PhC slab. [45] It turns out that this vertical confinement mechanism works much 
more effectively than the horizontal confinement mechanism by the PhC mirror after T > la, hence 
more radiation in the horizontal direction. For practical applications, however, we may redirect the 
horizontal propagating wave's energy into the vertical direction by employing grating couplers. [47] 

On the other hand, V tends to increase almost linearly with the slab thickness. It should be 
noted that V ~ l.^{X/nf at T ~ 720 nm, which is comparable with that of the widely used L3 
nanocavity (about 1.2 (A/n)'^). [11] We have also estimated the maximum achievable spontaneous 



UEM{t) 



(9) 



max{?7£;M(i)} 

Using Qtot (Eq. 2) and V (Eq. 9), Purcell's figure of merit can be written as [2] 




(10) 
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Fig. 6. (a) A PhC nanocavity is brought into contact with the underlying metal 
substrate. We assume realistic optical constants of gold at room-temperature, which 
is implemented using the single-pole Drude model in FDTD. (b) FDTD simulated 
electric- field intensity profiles (|Ep) of the dipole mode when the slab thickness is 
606 nm. Other structural parameters are as follows: Rm = 0.25a, R = 0.25a, Rp = 
0.05a, and a = 315 nm. (c) Optical properties of the dipole mode. Quality factor (Q), 
effective mode volume (V), and Purcell factor (Fp) derived from Q and V are plotted 
as a function of the slab thickness. Here, slightly different lattice constants (a) have 
been used for different slab thicknesses to keep the emission wavelength a constant 
at ~ 1.3;um. (d) The total electromagnetic energy contained in the cavity dissipates 
into two independent loss channels, one in the form of propagating radiation in air 
and the other in the form of absorption in metal. Here, radiation efficiency refers to 
the fraction of total dissipation into the radiation. 
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emission rate enhancement through the Pur cell factor Fp. [2] It is shown that, theoretically, Fp of 
over 100 can be achieved, a reasonably high value within the weak-coupling regime of the cavity 
QED. 

3.B. Threshold gain 

Qtot tells us how much threshold gain will be required to achieve lasing. According to the recent 
formulation by Chang and Chuang, [40,41] threshold gain, gtf^, is given by 

gth = 7^ : — , (lij 

i E ■ Vtot A 

where Ug^a is the material group index defined as d[ujna{uj)]/duj and na{u}) is the refractive index 
of the active region. F^; is the energy confinement factor defined as 

_ /y, c^^r Ui7;(r,t) +nA/(r,t) 

'- E — ~1- 7 — ^ 7 — r- v-*-^] 

Jy d^r UE{r,t) +UM[r,t) 

Here, the above volume integration for Va is taken over the volume of the active region and ue and 
UM are defined in Eqs. (5) and (6). Strictly speaking, in situations involving dispersive media, Te 
is a quantity that depends on time. However, in most cases, we can safely assume F^; as a constant 
for a given resonant mode. 

We take the PhC dipole mode as an example, whose electric-field intensity profiles are shown in 
Fig. 6(b). We choose the slab thickness to be 606 nm, roughly 2a. We find that Qtot and resonant 
wavelength are 1,032 and 1307 nm, respectively. Assuming the appropriate number of quantum 
wells as the gain material, F^; ~ 10% is not so difficult to achieve, from which g^^ is estimated to 
be ~ lOOOcm"^. This gth value is achievable by employing conventional InP/InGaAsP quantum 
wehs. [48] 

3.C. Radiation efficiency 

Because of the presence of the absorbing metal layer near the PhC cavity, energy contained in the 
cavity will be lost by two independent mechanisms: absorption loss in the metal (~ 1/Qabs) and 
radiation loss into air (~ l/Qrad)- Therefore, the total Q can be decomposed into the radiation Q 
and the absorption Q in the following manner. 

Ill 

^;tot ^;rad Vabs 

We define the radiation efficiency, 7?rad) as the ratio of total radiated power over the total dissi- 
pated power. This can be written in terms of quality factors in the following way. 

.rad-^ = l-^ (14) 

Here, Qtot can be easily estimated through field decay in the time-domain using Eq. (2). Estimating 
Qabs or Qrad requires additional volume integration or surface integration. For example, to calculate 
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absorbed power in the Drude metal, we should use the following volume integration. [42] 

i^abs(t) = I d^r weolm [e^(w)] (E(r, t) ■ E(r, t))^ (15) 
Jv 

Then, Qabs can be calculated by 

Qahs 

On the other hand, Qrad can be calculated by 

Qrnd ~ 

where -Prad(i) is given by 

Pradit) = j>^d\- (E(r, t) X H(r, t))^ (18) 

However, this type of surface integration usually requires a much larger computational domain 
since the surface of the integration should be located sufficiently far from the near-zone of the PhC 
nanocavity mode. Thus, we have adopted Eq. (16) for the calculation of the radiation efficiency. 

As expected, ?7i.ad is quite low (below 10 %) when T < 300 nm. However, we can bring up r^rad 
by simply increasing the slab thickness. It is found that r^rad becomes about 50 % at T = 600 nm. 
From this slab thickness, we can argue that the device's vertical radiation efficiency will begin to 
compete with the traditional PhC slab cavity suspended in air, by which at most 50 % of photons 
can be collected from the top (without the bottom mirror). In any case, such ?7rad values will be 
much higher than that of other classes of metallic nanocavities, where Qtot tends to be limited by 
Qabs. [34,36] 

3. D. Far- field emission 

Far-field directionality can be tuned by varying the slab thickness. Similar systematic optimization 
of the vertical collection efficiency has been reported for the nanowire cavity sitting on a flat metal 
surface. [49] For such a simple cavity geometry, a Fabry-Perot model can be used to optimize 
the far-field directionality. However, our PhC nanocavity involves complex geometrical features 
together with the zero group velocity dispersion along the z direction, making it difficult for us 
to develop a simple model as has been done for the PhC nanocavity suspended over the bottom 
mirror. Therefore, we have used 3-D FDTD and the near- to far-field transformation algorithm 
developed in our previous work. [15] Fig. 7 shows the evolution of far-field patterns as a function of 
the slab thickness. The result of T= 600 nm looks most promising; about 50 % of photon emitted to 
the top surface will be collected within ±30°. Fig. 7(b) shows that the vertical emission is linearly 
polarized along the x direction, whose direction of polarization has been determined by the two 
enlarged air-hole positions (See Figs. 2(b) and 6(b)). 

4. Degree of plasmonic eflfects: Plasmonicity 

As discussed in the previous section, the effect of the thick slab on the vertical confinement mech- 
anism can be understood from the viewpoint of a PhC fiber. On the other hand, the effect of the 
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Fig. 7. Here, we assume RT gold substrate and analyze the same dipole mode shown 
in Fig. 6(b). (a) Far-field emission profiles by varying the slab thickness from 500 
nm to 900 nm. (b) Polarization resolved far-field pattern for T = 600 nm. In the 
simulation, microscopic linear polarizers, one polarized along x direction and the 
other y direction, are assumed to scan over the hemispherical surface to measure 
\Ex\'^ and {Eyl"^, respectively. 
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bottom mirror may be understood by the method of image charges. [42] In the case of a PEC mirror, 
this scheme works perfectly — the fundamental resonant mode (even) in a PhC slab with thickness 
T is equivalent to the first-order mode (odd) in a PhC slab with thickness 2T. [6] However, in the 
case of a realistic metal mirror, the situation is not so simple. The existence of the evanescent field 
of the cavity mode in the metal does not preclude the possibility of interesting plasmonic effects 
within such a structure. 

So far, interesting aspects particular to SPP modes have been emphasized by many researchers. 
Two representative examples are the extremely large local electric-fields in a metal nano-antenna 
(hot spots) [50] and various meta-material engineering such as negative refractive index. [51] How- 
ever, we would like to emphasize that almost same performance results have been demonstrated 
from other engineered structures consisting purely of lossless dielectric media. [52, 53] Recently, 
Ishizaki and Noda have emphasized the similarity between the SPP wave and the surface state of 
a 3-D PhC possessing a 3-D PEG. [54] 

In fact, what makes metal metallic is the presence of negative permittivity rather than the metallic 
absorption originated from the imaginary part of the permittivity. One can easily show that any 
medium with e(a;) < cannot support any propagating planewave solution. [42] This result reminds 
us of light propagation behavior within the PBG. This observation suggests to us another viewpoint 
for metal — metal as a natural 3-D PBG material. In principle, we can carve and mold metal into 
an arbitrary geometrical shape. Interestingly, it has been known that even a tiny section (dimension 
< 10 nm) of metal does not loose its bulk optical properties, [55] which is in contrast to the case of 
an artificial 3-D PhC; at least several lattice periods are needed to function as a PBG material. [56] 

Another aspect of metal is that conduction electrons are strongly coupled to the electromagnetic 
fields. In fact, plasmons refer to the quantum mechanical eigenstates of coupled electron and photon 
states. [57] In the presence of damping (7^ > 0), the energy of plasmons will be converted into 
heat or mechanical energy. However, the damping mechanism is not an essential aspect of plasmons. 
Therefore, we can still think of the total energy of plasmons or the kinetic energy of conduction 
electrons. 

4. A. Definition of plasmonicity 

Here, we would like to find an expression for the kinetic energy density of electrons, Ukin(r, t), 
in metal. We denote the mass of an electron and the volume density of electrons as mg and N^,, 
respectively. Then, we can write the time-averaged kinetic energy of electrons per unit volume as 



■"kin(r, t) = -meNe{r, t) (v(r, t) ■ v(r, t))j, 



(19) 



The above expression can be rewritten in the following way. 




(20) 



(21) 
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In general, the current density vector J(r, t) is not a constant in time due to the presence of non- 
negligible damping. If we assume that damping is negligible within a few optical cycles, we can 
write J{r,t) w (1/2) ( J(r) exp(— iwt) + J*(r) exp(itjt) 



?^km(r) 



Nee^ 
nip 



J(r)|2 (22) 
a{u)\'\B{r)\' (23) 



Here, we have used Ohm's law J(r) = cj(a;)E(r). With the Drude model, we have 



7m - 

where to?, = NeC^ / {eonie)- Therefore, we obtain 



a(-) = (24) 



-Hn(r) = J-2^|E(r)P (25) 



im 

We can generalize the above result for the presence of damping such that 

,2 



nkin(r, t) = |-2^ (E(r, t) ■ E(r, t))^ (26) 

It is interesting to see that this kinetic energy density term can be found from ue and um (See 
Eqs. (7) and (8)). It should be noted that the role of kinetic energy density or kinetic inductance 
in developing a circuit model for plasmonics has been emphasized by Staffaroni, et al. [58] 

Now, we are in a good position to apply our understanding of Ukin to define a new quantitative 
measure of the degree of plasmonic character. Plasmonicity 11 can be defined the fraction of total 
kinetic energy over the total electromagnetic energy in the limit of zero damping as 

n ^ 2 X ""ZlT '"'''^ = 2 X , ^'^"'"^'l ^^ (m the limit 7m ^ 0) (27) 

Total EM energy Jy d^r {uE{r) + UM{r)) ^ im j \ j 

The reason we take the zero damping limit is that many actual experiments are performed at low 
temperature. [34, 35, 37] Even at room temperature, 7^ is smaller than w by at least one order of 
magnitude. Moreover, as we have emphasized in the above, the damping or the absorption is not an 
essential characteristic of metal. In this way, we can make the definition temperature-independent. 
The prefactor 2 is to scale up the maximum achievable 11 to 1, which will be shown in the following 
examples. 

In the following subsections, we shall examine 1) the SPP mode at a simple metal/dielectric 
interface and 2) the hybrid metal/PhC nanocavity, with which we shall show the validity of the 
definition of plasmonicity. 

4.B. The simple SPP mode 

We consider the SPP mode formed at a flat metal/dielectric interface. We assume two semi-infinite 
media that meet at the z =0 plane; a dielectric medium with ea for ^; > and the Drude metal 
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Fig. 8. We analyze the simplest surface plasmon polariton mode formed at a 
dielectric/metal interface, (a) Theoretical dispersion relation (u}-k) when = 
5.0 and e^o = 1-0, where we assmne the Drude model for the metal, €m{uj) = 
Coo — ^p/i^'^ + ijrn^)- (b) Plasmonicity is obtained by analytic calculation (solid 
curve) and compared with the result obtained using FDTD (square dots). 



with em(w) for z < 0. We also assume that the surface confined SPP mode propagates in the x 
direction. Then, the dispersion relation of the SPP mode is given by [59] 

kl=(-) (28) 
\cj Em + ea 

Figure. 8(a) shows the dispersion curve for = 5.0 and e^o = 1.0. The wavevector and the fre- 
quency are normalized by kp = Up/c and LOp, respectively. At small wavevectors, the curve increases 

1/2 

almost linearly; in fact, it asymptotically approaches the light line defined by = ck/ej . At large 
wavevectors, however, the curve rolls off and saturates at Wgp = i^p/(eoo + e^)^''^- Traditionally, the 
SPP mode at this large wavevector limit becomes the classical surface charge density wave and is 
considered as surface plasmon-like. In the opposite limit, of course, the mode is photon-like. Now 
let us quantify the degree of this plasmonic character. 
Using the definition of IT (Eq. (27)), we have 

"spp = ^Tjk; — 7^^^^^^ fr/tf — — ' (29) 



Jo Jo '^^ ^° \i) ^^"^ ' 

/ dx dz eoed {Ed ■ Ed)^ + / dx dz eofioo (E^ • Em)^ 
Jo Jo Jo Jo 



where E^ and E^ are electric-fields in metal and dielectrics, respectively. It is not so difficult to 
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Fig. 9. Here, we assume an ideal gold substrate by setting its damping constant to 
zero, (a) Quality factor of the same dipole mode as shown in Fig. 6. Since the ab- 
sorption in gold has been quenched, quality factors relate to radiation losses. Again, 
we tune the lattice constant (a) to approximately keep the emission wavelength at 
~ 1.3;um. (b) The degree of the plasmonic character, ^ plasmonicity" of the dipole 
mode as a function of the slab thickness. 



carry out the above integration analytically. Thus, we get 

nsPP = (30) 

One can show that IIspp reaches its maximum 1 in the limit of large wavevector where 
uj — LOp/ {eoo + ed)^^^ and — ^ —^d- The solid curve in Fig. 8(b) shows the analytical result from 
Eq. (30) . We have also performed FDTD simulations and the results are overlaid on the same graph, 
represented by square dots. 

4-C. Hybrid metal/PhC nanocavity 

Now, we return to the hybrid metal/PhC nanocavity and apply the definition of 11. It should 
be noted that we set ■jm = for all calculations in this subsection. First, we calculate Q of the 
same dipole mode as we vary the slab thickness. The lattice constant has been tuned to keep 
approximately the same emission wavelength of ~ 1.3/im. Figure 9(a) shows Q and the wavelength 
as a function of the slab thickness. Here, we would like to emphasize that Q obtained in this way 
is different from the radiation Q defined in Eq. (17), because the optical constants of gold have 
been changed. Though the change is small, it is large enough to produce a noticeable difference in 
Qrad because the reflectivity of the metal is not the same any more. The presence of a peak in Q at 
T = 300 nm seems to have the similar origin, the vertical resonance within the slab, as argued by 
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Tandaechanurat, et al. [43] Except for this feature, Q tends to increase almost linearly up to 2,000 
at T = 1,000 nm. 

We plot FDTD simulated Ilcav values in Fig. 9(b). When T > 600 nm, Ilcav is below 1%. Even 
when T = 275 nm, Ilcav is just above 4%. These values are much lower compared to those obtained 
by the SPP mode; only at very small wavevectors of k/kp < 0.2, such low values would be obtained. 
Therefore, our hybrid metal/PhC cavity mode is photonic-like rather than plasmon-like and the 
role of the metal mirror is not to generate plasmons but to serve as a mirror. 

5. Discussion and Conclusion 

We apply our definition of IT to the metal nanodisk cavity, which we have proposed recently. [38] 
Specifically, we take the two example modes, the monopole mode and the SPP-like mode, both 
emitting at ~ 664 nm (See Fig. 10). The SPP-like mode, as expected, shows a large 11 of 58.3% while 
the monopole mode shows 11 of 13.3 %. Therefore, we believe that the proposed formula provides 
reasonable quantitative measures for a wide range of metallic nanocavity designs. If metallic and 
dielectric structures of complex geometrical shapes play an important role in the formation of 
an electromagnetic mode of interest, a simple judgment based on the the mode profile can be 
misleading and does not provide any useful information on the internal properties of plasmons. 
However, the plasmonicty U requires simple volume integrations of energy densities. Thus, it can 
be easily applied to various situations regardless of how complex the cavity geometry is. 

There are two different routes of using metallic nanocavities for building coherent light emitters. 
One is to build an efficient light emitter, where the role of the metal is to serve as a mirror. 
One such example is our hybrid metal/PhC nanocavity. The other way is to build an efficient 
plasmonic source, also known as a SPASER (Surface Plasmon Amplification by Stimulated Emission 
of Radiation). [60] Such a (dark) plasmonic source would generate more plasmons rather than 
photons. In this context, it seems that the true SPASER should have very small radiation efficiency. 
However, as discussed in the previous sections, the conventional definition of radiation efficiency 
(Eq. (14)) could be misleading since it contains the effect of the irreversible damping process and the 
metallic absorption is not an essential property of plasmons. The concept of plasmonicity defined 
here may provide a useful guideline for the optimization of the SPASER source. By the same token, 
it could serve equally well for the development of efficient photon sources based on a metallic 
mirror. In both applications, however, increasing Qtot will be critical to achieve the condition for 
the stimulated emission. 

Acknowledgments 

The authors would like to acknowledge support from the Defense Advanced Research Projects 
Agency under the Nanoscale Architecture for Coherent Hyperoptical Sources program under grant 
#W911NF-07-l-0277 and from the National Science Foundation through NSF CIAN ERG under 
grant #EEC-0812072. 



22 



{a) SPP-like mode 





SPP-like 


Monopole 


^(nm) 


664 


663 


Q 


317000 


1313 


n 


0.583 


0.133 




Fig. 10. We take an example from our previous paper, a metal-clad nanodisk. [38] 
The diameter of the dielectric disk is 220 nm. Here, two representative modes are 
analyzed, (a) one is a SPP-like surface confined mode and (b) the other is a photonic- 
like monopole mode. Their emission wavelengths, quality factors, and plasmonicities 
are obtained through FDTD simulations. All values correspond to ideal silver with 
zero damping. 
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